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The thermodynamics of a first-order chiral phase transition is considered in the presence 
of spinodal phase separation using the Nambu-Jona-Lasinio model in the mean field 
approximation. We focus on the behavior of conserved charge fluctuations. We show 
that in non-equilibrium the specific heat and charge susceptibilities diverge as the system 
crosses the isothermal spinodal lines. 



1. Introduction 

The search for the critical end poin t (C EP) is one of the central issues in strongly 
interacting hot/dense QCD matter Ui2i. It is of particular interest to identify the 
position of the CEP in the phase diagram and to study general properties of ther- 
modynamic quantities in its vicinity. Modifications in the magnitude of fluctuations 
or the corresponding susceptibilities can be considere d as a possible signal for de- 
confinement and chiral symmetry restoration | 3 | 4 | 5 | 2 E)| j n ^jg con t e xt, fluctuations 
related to conserved charges play an important role since they are directly accessible 
in experiments^. 

The enhancement of the baryon number fluctuations could be a clear indication 
for the existence of the CEP in the QCD phase diagram. However, finite density 
fluctuations along the first-order transition appear under the assumption that this 
transition happens in equilibrium. This is modified when there is a deviation from 
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equilibrium In this contribution we briefly show that the enhanced baryon number 
density fluctuations is a signal for the first-order phase transition in the presence of 
spinodal decomposition. 

2. Quark number susceptibility and spinodal instabilities 

In a non-equilibrium system, a first-order phase transition is intimately linked with 
the existence of a convex anomaly in the thermodynamic pressure There is an 
interval of energy density or baryon number density where the derivative of the 
pressure, dP/dV > 0, is positive. This anomalous behavior characterizes a region 
of instability in the (T, n g )-plane which is bounded by the spinodal lines, where 
the pressure derivative with respect to volume vanishes. The derivative taken at 
constant temperature and that taken at constant entropy, 

dP\ fdP\ 

^j) -0 and — =0, (1) 



dVJ T \9V /S 

define the isothermal and isentropic spinodal lines respectively. 

If the first-order phase transition takes place in equilibrium, there is a coexistence 
region, which ends at the CEP. However, in a non-equilibrium first-order phase 
transition, the system supercools/superheats and, if driven sufficiently far from 
equilibrium, it becomes unstable due to the convex anomaly in the thermodynamic 
pressure. In other words, in the coexistence region there is a range of densities and 
temperatures, bounded by the spinodal lines, where the spatially uniform system is 
mechanically unstable. Spinodal decomposition is thought to play a dominant role in 
the dynamics of low energ y nuclear collisions in the regime of the first-order nuclear 
liquid-gas transition^ ' '. Furthermore, the consequences of spinodal decomposition 
in the connection with the chi ral and deconfinement phase transitions in heavy ion 
collisions have been discussed I 9 | 1U | H | 

In Fig. [ljleft we show the evolution of the net quark number fluctuations along 
a path of fixed T = 50 MeV calculated in the Nambu-Jona-Lasinio (NJL) model 

iini 

in the mean field approximation . When entering the coexistence region, there 
is a singularity in \q that appears when crossing the isothermal spinodal lines and 
where the fluctuations diverge and the susceptibility changes sign. In between the 
spinodal lines, the susceptibility is negative. Consequently, this implies instabilities 
in the baryon number fluctuations when crossing from a meta-stable to an unstable 
phase. The above behavior of \q IS a direct consequence of the thermodynamics 
relation 

-) =- nl ~- (2) 
dV) T V Xq [> 

Along the isothermal spinodals the pressure derivative in Eq. @ vanishes. Thus, 
for non- vanishing density n q , \q must diverge to satisfy Furthermore, since the 
pressure derivative 8P/dV\T changes sign when crossing the spinodal line, there 
must be a corresponding sign change in x g , as seen in Fig. [TJleft. 
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Fig. 1. (Left) The net quark number susceptibility at T = 50 MeV as a function of the quark 
number density across the first-order phase transition. (Right) The net quark number susceptibility 
in the stable and mcta-stablc regions. 



In Fig.QJright we show the evolution of the singularity at the spinodal lines when 
approaching the CEP. The critical exponent at the isothermal spinodal line is found 
to be 7 = 1/2, with Xq ~ (a* - aO~ 7 , whilc 7 = 2 /3 at the CEP®. Thus, the singu- 
larities at the two spinodal lines conspire to yield a somewhat stronger divergence 
as they join at the CEP. The c ritical region of enhanced susceptibility around the 
TCP/CEP is fairly small QMS while in the more realistic non-equilibrium system 
one expects fluctuations in a larger region of the phase diagram, i.e., over a broader 
range of beam energies, due to the spinodal instabilities. 

The rate of change in entropy with respect to temperature at constant pressure 
gives the specific heat expressed as 



= TV 



2s 

XTT XfJ.T 

n q 



(3) 



The entropy xtt and mixed x^t susceptibilities exhibit the same behaviors as 
that of Xq shown in Fig. QJleft. Thus Cp also diverges on the isothermal spinodal 
lines and becomes negative in the mixed phase Q. It was argued that in low energy 
nuclear collisions the negative specific heat could be a signal of the liquid-gas phase 
transition Its occurrence has recently been reported as the first experimental 



evidence for such an anomalous behavior 



3. Conclusions 

We have shown that in the presence of spinodal instabilities the net quark number 
fluctuations diverge at the isothermal spinodal lines of the first-order chiral phase 
transition. As the system crosses this line, it becomes unstable with respect to spin- 
odal decomposition. The unstable region is in principle reachable in non-equilibrium 

a The specific heat with constant volume, on the other hand, continuously changes with n q and 
has no singularities on the mean-field level. 
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systems that is most likely created in heavy ion collisions. Consequently, large fluc- 
tuations of baryon and electric charge densities are expected not only at the CEP 
but also when system crosses a non-equilibrium first-order transition. 
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